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STRESS  AND  THERMAL  STRESS  COMPENSATION  IN  QUARTZ  SAW  DEVICES 


GKNERAL  PRESENTATION 


Surface  acoustic  waves  are  sensitive  to  different  physical  quantities  such  as 
temperature,  forces,  accelerations, ...  [1-8].  The  external  perturbations  induce  static 
or  quasi-static  stresses  and  strains  in  the  crystal,  which  by  non-linear  coupling 
modify  the  wave  velocity  and  therefore  induce  frequency  shifts  in  a  SAW  device.  It 
is  necessary  to  minimize  these  effects  in  the  case  of  SAW  oscillators  since  they  are  at 
the  origin  of  most  of  frequency  instabilities.  (On  the  other  hand,  high  sensitivities 
can  be  used  also  for  sensor  applications). 

Static  temperature  effects  may  be  avoided  by  the  choice  of  a  temperature- 
compensated  cut  like  the  classical  ST-cut. ' 

The  effects  of  forces  depending  closely  on  the  geometry  of  the  device  (shape  and 
fixation  conditions),  one  may  think  to  study  this  problem  for  each  particular  device. 
Another  way  which  is  followed  in  the  present  work  is  to  investigate  the  possibility  of 
stress  compensation  by  using  the  crystal  anisotropy  properties.  (Fig.l) 


COMPENSATION  OF  STRESS  EFFECTS 
IN  SAW-DEVICES 


Conventional  cut 
(e.g:ST,X) 


Optimized  plate  shape 
and 

fixation  conditions 


Using  crystal  anisotropy 

to  find 

intrinsic  compensation 

to  certain  kinds 

of  stresses 

NEW  CUTS 

Fig.l :  The  problem  of  stress  sensitivity  in  SAW  devices 
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In  this  approach  [9,10]  the  effects  of  stresses  are  considered  as  the  tensorial 
equivalent  of  temperature  effects,  and  under  some  hypothesis  stress  sensitivity 
coefficients  can  be  defined  'ivhatever  the  origin  (thermo-elastic  effects,  external 
forces)  might  be.  Therefore  a  general  relation  between  the  wave  velocity  and  the  six 
independent  components  of  the  quasi-static  stress  tensor  can  be  given. 

Stress  compensation  will  be  obtained  if  the  contributions  of  the  main  stress 
components  cancel  each  other,  and  this  may  happen  for  a  given  distribution  of 
stresses  when  a  suitable  SAW  quartz  cut  is  found  minimizing  the  velocity  shifts. 
However  as  a  compensaj'on  seems  unrealistic  in  the  most  general  case,  the  problem 
will  be  treated  first  in  the  case  of  planar  stresses  in  thin  plates,  which  reduces  the 
number  of  stress  components  to  three.  In  addition,  a  symmetrical  stress  distribution 
will  be  considered,  corresponding  to  equal  axial  stresses  and  negligible  shear 
component.  In  that  configuration,  stress  compensation  means  that  the  two  in-plane 
axial  stresses  induce  equal  contributions  with  opposite  signs.  This  approach  was 
used  by  Eernisse  for  bulk  acoustic  wave  resonators  [11]  and  led  to  the  development 
of  the  famous  SC-cut. 

The  previous  analysis,  which  only  takes  into  account  the  contribution  of  in-plane 
stresses,  has  been  extended  using  the  same  formalism  to  the  investigation  of 
flexural  sensitivity  of  quartz  cuts.  In  the  case  of  pure  bending,  the  SAW  (which  is 
located  close  to  the  surface)  will  be  submitted  to  compressive  in-plane  stresses  only  ; 
a  symmetrical  distribution  of  such  bending  forces  will  therefore  be  equivalent  to  a 
symmetrical  distribution  of  compressive,  in-plane  stresses  and  will  not  induce  any 
frequency  shift  in  a  SAW  cut  compensated  for  planar  isotropic  stresses. 

The  present  work  shows  that 

-  several  quartz  cuts  exist,  exhibiting  both  zero  first  order  temperature  coefficient 
and  zero  stress  coefficient  for  planar  isotropic  stresses,  according  to  the  model. 

-  among  these  cuts  exist  two  singly-rotated  quartz  plates,  and  experimental  tests 
were  performed  to  check  temperature  and  stress  sensitivity  of  these  particular  cuts  ; 
results  are  compared  to  theoretical  predictions  obtained  by  analytical  methods  and 
finite  element  analysis. 

-  theoretical  calculations  on  flexural  sensitivity  of  the  new  cuts  (SAW  propagating 
on  a  bent  cantilever  beam,  on  a  plate  submitted  to  normal  accelerations  or  on  a  thin 
rectangular  anisotropic  plate)  predict  a  possibility  of  bending  forces  compensation  if 
a  symmetrical  distribution  of  bending  stresses  is  applied  to  the  device. 


5 


1. 


DKRIVATION  OF  STRESS  SiilNSlTl VITY  COEFFICIENTS  FOR 
SURFACE  ACOUSTIC  WAVES 


a)  Wave  propagation  in  a  prestrained  medium 

The  calculation  is  based  on  the  model  proposed  by  Tiersten  [12,  13]  for  wave 
propagation  in  a  prestrained  medium. 


The  general  equation  of  a  finite  amplitude  wave  of  displacement  Uj  propagating  on  a 
medium  with  specific  mass  po  in  the  natural  state  submitted  to  a  mechanical  bias 
[12]  can  be  written  as 


p  u  = 

^0  I 


a  - 

—  {A.  . 
da 


du 

r 

da 


(la) 


_  da^ 

— 1  ~  0  on  a  free  surface 
s 


(lb) 


where  S/dak  refers  to  the  natural-state  coordinates,  Nk  is  the  cosine  director  of  the 
normal  to  the  surface.  Aiskr  are  the  modified  elastic  coefficients  resulting  from  the 
nonlinear  coupling  between  the  static  stresses/strains  and  the  high  frequency  wave. 
These  coefficients  can  be  written  in  the  form 


^iskr  ^iskr 


// 


tskr 


(2) 


where  Ciskr  are  the  regular  second  order  elastic  constants,  and  Hiskr  appears  as  a 
perturbation  due  to  the  bias,  given  by  the  relation 


(3) 


5ik  is  the  Kronecker  symbol,  Ciskruv  are  the  third  order  nonlinear  elastic  constants, 
Tsr,  Suv  and  317, •  /  dae  represent  respectively  the  static  stresses,  strains  and 
displacement  gradients  characterizing  the  mechanical  bias.  (Summation  convention 
on  repeated  indices  is  assumed). 


b)  General  expression  of  SAW  velocity 

The  mechanical  displacement  corresponding  to  a  SAW  propagating  in  the  ai 
direction  with  a  velocity  Vo  at  an  angular  frequency  coo  on  an  anisotropic  medium  in 
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the  absence  of  perturbation,  taking  piezoelectricity  into  account,  is  expressed  by 


u 


0 

I 


q(r)a^ 

- 

•*  0  V 

b  (r)  A  e  “ 

^  I  r 
r  =  1 


(j 


e 


(4) 


The  reference  axes  (aj)  of  the  natural  state  are  presented  in  Fig.  2  with  respect  to  the 
crystal  plate. 


Fig.  2  :  Reference  axes  for  a  SAW  propagating  on  a  semi-infinite  crystal  plate 

Arbiir)  represents  the  amplitudes  of  the  4  partial  waves  (r  =  1-4)  of  the  unperturbed 
SAW.  q(r)  is  the  depth  penetration  coefficient  of  the  r-th  partial  wave  in  the  a2- 
direction. 


Eq.  (4)  may  be  expressed  in  a  closer  form  by 


u 


0 

I 


4 


1(0  ({  -  («.  a,)/V 

0  C  C  (J 

e 


where 


n^{r)  =  1  Ai,(r)  =  q{r)  n^{r)  =  0 


(5) 


Applying  a  perturbation  method  [13, 14]  to  the  SAW  propagating  on  the  constrained 
crystal  (of  volume  V)  the  relative  velocity  shifts  A V/Vo  can  be  expressed  as 


A’/ 


V 

0 


au?  dal 

-i-irr  —  d» 

da  da 


2p  (0^ 
‘^11  o 


o  o 
u.  u 

I  I 
■J 


dv 


(6) 
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In  this  approach,  piezoelectricity  is  taken  into  account  only  to  determine  the 
unperturbed  SAW  expressed  as  a  sum  of  4  partial  waves.  The  perturbation  term 
Hiskr  is  the  same  as  in  the  non-piezoelectric  approximation. 


In  the  most  general  case,  the  perturbation  terms  Hiskr  depend  on  the  space 
coordinates  if  the  distribution  of  the  static  stresses  and  strains  is  inhomogeneous, 
and  a  particular  distribution  of  stresses  is  associated  wi  th  each  particular  boundary 
condition.  If  considering  the  case  of  an  homogeneous  stress  and  strain  distribution 
eq.  (6)  takes  the  simpler  form : 


AV  - 

Y  iskr  iskr 

0 

where  Uiskr  depends  only  on  the  unperturbed  SAW  parameters 


U 


iskr 


1 

9 

2p  0}" 
^0  0 


t  au?  au° 

—  —  du 

Jm  aa  Ba 
s  r 


(7) 


(8) 


Assuming  a  plate  thickness  much  greater  than  the  penetration  depth  of  the  SAW, 
the  integration  of  eq.  (8)  over  02  0]  yields  a  similar  expression  as  in  [14],  except 

that,  taking  piezoelectricity  into  account,  summation  is  made  over  4  partial  waves 

4 

A  b,(p)  A  b  (q)  n  ip)  n  {q)/{q  -q) 

y  ^  ^  M=1 _ 

2p  4  A  bAp)  A*  bAq)* 

''00  -y  P  J  H  J  ^ 

p.q=i  (Ip  - 


c)  Stress  sensitivity  coefficients 


The  perturbating  coefficients  Hiskr  implicitly  contain  the  dependence  between  the 
velocity  and  the  static  bias  (eq.  3)  in  term  of  stresses,  strains,  and  mechanical 
displacements.  The  strains  can  be  expressed  as  function  of  the  stresses  by  means  of 
the  compliance  elastic  constants  Suvpq- 
Ciskruv  Suv  is  written  : 


C 


iskr 


s 


uu 


=  c 


iskruo 


S 

uupq 


(10) 


The  two  last  terms  of  eq.  (3),  bUi!  da(  and  bUk!  bat,  cannot  be,  in, general,  directly 
expressed  as  a  function  of  the  strains  Suy. 

But  for  solving  the  problem  of  the  static  deformations,  it  is  necessary  to  apply  the 
fixation  conditions  of  the  plate.  A  condition  of  a  fixed  point  at  the  origin  (ai  =  02  = 
as  =  0)  is  first  given.  A  condition  of  no  rigid  rotations  of  the  crystal  plate  around  the 
fixed  point  may  be  introduced  without  loss  of  generality.  This  second  condition 
implies 


at/  at/, 

I  ^  £ 

ac^  "  aa, 


for  Oj  =  a,  =  =  0 


with  as  consequence 


—  =  S  for  a  =  a  =  a  =  0 


Assuming  that  the  SAW  propagates  in  a  small  volume  close  to  the  origin  eq.  (12)  can 
be  considered  to  be  satisfied  in  this  part  of  the  crystal  plate.  Expanding  Sit  in  terms 
of  the  stresses  Tpq,  eq.  (12)  may  be  expressed  close  to  the  origin  as 


—  =  s.,  T 
da^  w 


Combining  eqs.  (10)  and  (13)  into  eq.  (3),  the  velocity  shift  ^VIVq  is  obtained  in  the 
final  form 

AV  - 

—  =  T  (8,6  8  +  C.  ,  s  +  C,  ,  s.,  +  C.  ,  }  t/.  ,  (14) 

Y  pq  ^  ik  sp  rq  iskruo  uvpq  Cskr  itpq  isCr  k€pq*  iskr 

0 


or  in  a  condensed  form 


—  =  r  ‘q 

V  M 
0 


where  ^^apq  represent  the  different  stress  coefficients.  They  have  the  same  symmetry 
as  Tpq,  i.e.  there  are  six  independent  stress-sensitivity  coefficients. 


d)  Numerical  example  for  quartz 


As  it  may  be  seen  in  eq.  (14),  the  stress  sensitivity  coefficients  depend  only  on  the 
crystal  anisotropy,  and  are  independent  of  the  plate  geometry. 
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They  are  computed  by  tensorial  contraction  on  two  independent  terms  : 

-  the  first  term  is  a  contraction  on  2nd  order  and  3rd  order  elastic  constants  : 

CijkC>  Ciskruv 

-  the  second  term  Uiskr  is  a  combination  of  terms  depending  only  on  the 
unperturbed  SAW  parameters. 

First,  the  elastic  constants  are  "rotated"  and  expressed  in  the  set  of  axes  parallel  to 
the  surface,  as  defined  on  Figs.  2  and  3.  Then  the  surface  wave  parameters  are 
computed,  taking  piezoelectricity  into  account,  using  a  very  classical  algorithm, 
close  to  Campbell  &  Jones  method  [15].  This  gives  the  term  Uiskr,  which  is 
contracted  with  the  first  bracket  of  eq.  (14)  to  give  the  stress  sensitivity  coefficients. 

An  example  of  numerical  results  for  quartz  crystal  is  given  on  Fig.  4  where  the  six 
independent  stress  coefficients  are  plotted  on  a  contour-line  map,  as  a  function  of  cut 
angles  0,  and  propagation  direction  of  the  SAW  These  angles,  defined 
according  to  IEEE  standard  [16]  are  represented  on  Fig.  3,  This  corresponds  to  a 
(YXwf t)4>0^  plate,  SAW  propagation  in  the  f-direction. 


Fig.  3  :  Definition  of  the  cut  angles  0, 0  and  propagation  direction  'T 


IQ 

Correspondance  with  Euler  angles  [17]  is  the  following : 
A-Euler  = 

PEuler  = 
eEuler  =  'I'ikkk 


Examples 

(for  a  detailed  discussion  of  the  definition  of  cut  angles  and  symmetry  properties  of 
quartz,  see  Appendix  1). 


cut 

Euler 

IEEE 

Y-cut,  X  prop. 

(.\  =  0,  p  =  90°.  0  =  0) 

(cj)  =  0, 0  =  0,  'F  =  0) 
equivalent  for  SAWs  to 
((J)  =  0,0  =  180°,  'F  =  0) 

ST-cut,  X  prop. 

(X  =  O.p  =  132,75°. 
0=0) 

(4,  =  0.  0  =  42,75,  ^  =0) 
equivalent  for  SAWs  to 
(4  =  0,0  =  222,75,  ’F  -  0) 

X-33  cut 

(X  =  90°.  p  =  90°. 

0  =  33,44°) 

(4  =  90°,  0  =  180°,  T  =  33,44°) 
equivalent  for  SAWs  on  quartz  to 
(4  =  -90°,  0  =  0,  'F  =  33,44°) 

Numerical  computation  is  performed  using  FORTRAN-77  programs  on  a  DATA 
GENERAL  MV  1500/100  computer.  Elastic,  piezoelectric  and  dielectric  constants 
are  taken  from  Slobodnick  [17],  Bechmann  [18]  and  Thurston  [19]. 

A  complete  mapping  (Fig.  4)  0°  <  0  <  90°  (5°  steps)  and  -60°  <  <  60°  (5°  steps) 

requires  19  X  25  values  and  takes  about  15  minutes  of  CPU  Time. 

Values  of  cut  angles  non  represented  on  Fig.  4  may  be  found  according  to  symmetry 
properties  of  SAWs  on  quartz  (see  Appendix  1) 

(4,0, -40  =  ‘u^(-4,8.¥) 

.-*1 

‘a  (4, -0,40  =  *a  {ldO'’-4,0,4>) 

p(f^  pq  -r  »  » 


(17) 


a{2.2 


8(>-50-4(>-30-2(>-l0  0  10  20  30  40  50  60 


-10 
2 

3 

4 

5 

6 
7 
0 

^-60-50-4 


10  0  10  20  30  40  50  80 


a  2.3 


L-i-e-. 


60-5(>-40-30-2(>-10  0  10  20  30  40  50  80  ^  -6(>-5(>-4O-3O-2()-I0  0  r0'20  30  40  50  60 


vs.  s.  * 


Dill 


-6(>-5(>-40-30-20-10  0  10  20  30  40  50  60  ^  -60'50-40-30-20-10  0  10  20  30  40  50  60 


Fig.  4 

Example  of  numerical  results  for  stress  coefficients 
Contour-line  mapping  for  SAWs  on  Quartz,  propagation  angle  ip = 30® 
Piezoelectricity  of  the  crystal  is  taken  into  account 
units :  10'*^  Contour  lines  interval :  5  X  10’*^  m^.N'‘ 

solid  lines - sa  >  0  dotted  lines ....  sq  =  0  broken  lines  —  <  0 
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2.  QUARTZ  CUTS  COMPENSATED  FOR  BOTH  1st  ORDER  TEMPE¬ 
RATURE  EFFECTS  AND  PLANAR  ISOTROPIC  STRESSES  EFFECTS 


In  this  section,  the  model  of  stress  sensitivity  coefficients  is  applied  to  the  problem  of 
planar,  isotropic  stress  compensation.  Sensitivity  to  static  temperature  effects  has 
also  to  be  taken  into  account,  to  define  SAW  quartz  cuts  with  low  stress  and 
temperature  sensitivity. 

a)  Symmetrical  in-plane  stress  distribution  in  a  thin  plate 

The  case  of  a  thin  plate  of  quartz  with  the  main  surfaces  free  of  stresses  will  be 
considered  now.  In  this  case,  there  is  a  planar  stress  distribution  of  components  Tn, 
T33  and  T31  expressed  in  the  set  of  axes  parallel  to  the  surface  (Figs.  2  and  3). 

An  additional  synunetry  condition  will  be  assumed  :  Tn  =  T33,  T13  =  0.  This 
corresponds  to  a  planar  isotropic  stress  distribution  which  for  instance  appears  in 
the  problem  of  thermal  stresses  in  a  circular  plate.  With  these  assumptions,  the 
velocity  shifts  depend  only  on  a  single  stress  coefficient  sa 

AV  .  — 

—  =  "a.T  (18) 

0 

where  T  =  Tn  =  T33andsa  =  San  +  ®a33 

Compensation  of  planar  isotropic  stresses  using  quartz  anisotropy  consists  in 
finding  SAW  quartz  cuts  (4),9,qf)  satisfying 

(\  +  ‘033)  (4>,d.w)  =  0 


b)  1st  order  static  temperature  effects 


Static  temperature  sensitivity  will  be  defined  here  using  the  temperature 
coefficients  of  fundamental  elastic  constants  according  to  H.F.  Tiersten  and  B.K. 


1,12]. 


he  same  formalism  as  in  sec.  1,  the  first  order  temperature  coefficient  0a  is 
for  a  pure  thermal  bias  in  a  free-expanding  plate  as 


V 


=  (0  -  G  ) 

O 


(20) 


where 
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'^Q  =  {c  ,  e  +  c  .e.  +  c-  t  , )}  u. , 

likruu  uu  qskr  iq  nqr  kq  ^0  iskr  ‘  tskr 


and  (B-d,))  is  the  temperature  elevation  with  respect  to  the  reference  state. 


In  eq.  (21)  are  the  coefficients  of  linear  expansion,  d(CiskrVciQ  are  the  first 
derivatives  of  the  fundamental  elastic  constants,  and  (7iskr  is  the  combination  of 
SAW  parameters  defined  in  the  perturbation  method  for  an  homogeneous 
distribution  of  stresses  [10]. 


The  loci  of  stress-compensated  cuts  are  defined  according  to  eq.  (19).  Taking  the 
propagation  angle  ip  as  a  parameter,  these  loci  are  represented  on  Fig.  5  in  the  (0,$) 
\  plane  in  solid  lines  for  different  propagation  directions  :  ip  =  30“  (Fig.  5a),  vp  =  45° 
(Fig.  5b).  On  the  same  graphs  are  represented  in  broken  lines  the  loci  of  zero  first- 
order  temperature  coefficients. 


^  =30®  ■  qj  =45® 


Solid  lines :  locus  of  quartz  cuts  compensated  for  planar  isotropic  stresses  effects 
dotted  lines :  locus  of  zero  first  order  temperature  coefficient 
cut  angles  0,  $  and  propagation  direction  \p  correspond  to  a  (YXw€t)<l)0vp  plate 
SAW  propagation  in  the  f-direction,  according  to  IEEE  standard  [16] 
Piezoelectricity  taken  into  account 
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The  intersections  between  the  two  loci  define  quartz  cuts  satisfying  both  conditions 
of  stress  and  temperature  compensation.  Such  intersections  have  been  found  to  exist 
for  angles  between  20°  and  60°  approximately. 

In  the  case  of  surface  waves,  the  additional  parameter  qr  (if  compared  to  bulk  waves) 
allows  a  degree  of  freedom  in  the  choice  of  stress-  and  temperature-compensated 
cuts.  As  a  first  criterium  to  choose  between  the  different  intersection  points  shown 
on  Fig.  5,  it  is  interesting  to  look  for  singly  rotated  quartz  plates  (<1>  =  0)  since 
machining  tolerances  are  easier  to  obtain  than  for  doubly  rotated  plates. 

The  loci  of  temperature-  and  stress-compensated  cuts  are  represented  on  Fig.  6  for 
singly  rotated  cuts  (4>  =  0). 


0  Singly  rotated  Quartz  cuts  ^=^0 


-  zero  planar  isotropic  stress  coefficient 

Fig.  6 

Solid  lines :  locus  of  quartz  cuts  compensated  for  planar  isotropic  stresses  effects. 

Broken  lines :  locus  of  zero  first  order  temperature  coefficient 

Cut  angles  <1),6  and  propagation  direction  tp  correspond  to  a  (YXw^t)  4)6^  plate 
SAW  propagation  in  the  ^-direction  according  to  IEEE  standard  1949  [16] 
Piezoelectricity  taken  into  account 
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Two  intersections  exist 

-  for  ly  close  to  30®  (point  PSd) 

-  for  ip  close  to  45®  (point  P45) 

This  last  cut  (P45 :  (|>  =  0, 9  =  45°,  ip  =  45°)  is  close  to  B.K.  Sinha's  STC-cut  [20] 

(<^  =  0, 9  =  41.8°,  ip  =  46.89°). 

In  order  to  check  the  model  presented  here,  experimental  work  has  been  performed 
on  these  two  cuts,  as  exposed  now. 


3.  EXPERIMENTAL  RESULTS 


a)  Experimental  method 

Experiments  are  performed  on  quartz  delay  lines  (100  and  150  finger  pairs, 
transducer  period  :  34  pm,  time  delay  =  1.5  ps)  built  on  circular  plates  (diameter  : 
22  mm,  thickness  :  1  mm  or  2  mm).  The  device  is  used  in  an  oscillator  operating  at 
100  Mhz. 

Temperature  compensation  is  checked  by  plotting  the  frequency  vs  temperature 
curve  between  -20°C  and  +50°C.  The  propagation  angle  ip  is  adjusted  to  keep  the 
inversion  point  of  the  frequency/temperature  curve  between  18°C  and  25°C. 

When  temperature  compensation  has  been  found,  the  delay  line  is  submitted  to 
mechanical  stresses  with  an  apparatus  represented  on  Fig.  7. 

The  circular  plate  is  submitted  to  a  single  or  double  diametral  compression  obtained 
by  two  pairs  of  knife-edges  acting  at  90°  to  each  other. 

The  plate  may  be  rotated  around  its  axis,  and  the  force-frequency  effect  is  measured 
for  each  azimuthal  angle  P  between  the  propagation  direction  of  the  SAW  and  the 
direction  of  the  force  (Fig.  7). 

The  experimental  procedure  is  as  follows ; 

The  knife-edges  are  removed,  and  the  plate  is  held  by  a  clamp  until  a  new  angular 
position  is  set.  The  plate  is  then  submitted  to  a  diametral  pre-compression  of  about 
IN.  The  clamp  is  removed,  so  that  the  plate  is  submitted  only  to  the  compressive 
forces.  An  additional  force  between  IN  and  ION  is  set,  and  the  frequency  shifts  are 
measured  with  respect  to  the  pre-strained  state.  The  azimutal  angle  is  measured 
without  contact  by  optical  means. 
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The  force-frequency  effect  in  a  circular  disk  may  be  derived  from  the  stress  sensivity 
coefficients  and  from  the  classical  theory  of  elasticity  [21]. 


Fig;7  :  Experimental  setup  for  single  and  symmetrical  diametral  compression  of  a 

circular  quartz  delay  line 


The  stress  components  at  the  center  of  an  isotropic  disk  (diameter  d,  thickness  t) 
submitted  to  a  single  diametral  compression  by  a  force  P  are  expressed  as 


=  -  —  ,  T" .,  =  +  —  ,  Tl,  =  0 
ndt  '1^  ndt  1*1 


(23) 


where  upper  index  "o"  refers  to  a  set  of  axes  (aiH.aaH)  parallel  to  the  direction  of  the 
force  (Fig.  7). 


After  rotation  to  a  set  of  axes  (ai.as)  parallel  to  the  direction  of  the  SAW,  the 
relative  velocity  shifts  may  be  expressed  in  the  final  form 


AF  2P  /  .  ,  ,  \ 

F  "  "  iidi  (^‘"‘ll  ‘“33^  ^^‘‘'ll  ~  '"33^  ~  ^  °13 


Eq.  (23)  shows  that  the  angular  dependence  of  the  frequency  shifts  is  sinusoidal, 
with  an  offset  value  of  (-2P/ndt)  (Son  +  8033)  corresponding  to  the  case  of  a 
symmetrical  compression  with  two  orthogonal  forces  of  magnitude  P/2. 


b)  Results  for  the  cuts  ((1)  =  0,  -13°  <  9  <  -2.5°,  29°  <  u/  <  33°) 
and  (d)  =  0.  37°  <  9  <  47°.  42°  <  w  <  50°). 

Temperature  effects 

Singly-rotated  cuts  in  the  vicinity  of  the  P30  cut  were  tested  with  cut  angles 
9  =  -2.5°,  -5°,  -7°,  -10°,  -13°,  and  the  propagation  angles  qf  corresponding  to  an 
inversion  point  at  room  temperature  were  found  at  32.5°,  31.5°,  31°,  30°,  29° 
respectively  (Figs.  8  and  9).  As  a  reference,  the  theoretical  curve  (-34x10-9  ppm.K-2) 
for  the  classical  ST,X  cut  is  shown.  It  may  be  seen  that  the  thermal  behavior  of  the 
new  cuts  is  very  similar  to  the  classical  ST,X  cut.  In  the  vicinity  of  the  P45  cut,  four 
singly-rotated  cuts  (4>  =  39°,  qf  =  46° ;  6  =  33°,  ip  =  44.3°  ;  9  =  41.8°,  ip  =  48,2°  ; 
9  =  45°,  ip  =  47°)  were  tested.  Results  appear  on  Figs.  10  and  11. 

Temperature-compensated  cuts  were  experimentally  found  close  to  the  theoretical 
values  (jtptheory "  Wexpl)  <  0.5°)  in  the  vicinity  of  the  P30  cut.  The  difference  between 
theory  and  experiment  is  greater  (about  2°)  in  the  vicinity  of  the  P45  cut,  and  the 
corrections  introduced  by  the  piezoelectric  model  do  not  cancel  this  discrepancy 
existing  also  in  the  non-piezoelectric  approximation  used  before. 
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EXPERIMENTAL  FREQUENCY-TEMPERATURE  CURVES 


Fig.  8  :  Experimental  static  temperature  characteristics  for  the  cuts 
(4)  =  0,  -13“  <  0  <  -2.5“,  29“  <  ip  <  33“) 

Comparison  with  the  classical  ST,X  cut 


4^=0 


I  I - ^ — m - : — - - 1 

as  as  30  3a  34 


-  STRESS  0OEFRCIENT  (unils  mVN"') 

.  TEMPERATURE  COEFRCtENT  (units  :ppm.K'’) 

Fig.  9 : 

Theoretical  contour  line  maps  of  temperature  and  planar  isotropic  stress  coefficients 
Thess  curves  .are  used  as  a  guide  line  to  experiment  real  quartz  cuts 
in  the  vicinity  of  the  intersection  of  both  zero-stress  and  zero-temperature  effects 

cut  angles :  <^  =  0,  -13“  <  0  <  -2.5“ 
propagation  direction :  26“  <  ip  <  34“ 

«  =  real  cuts,  experimentally  found  temperature-compensated 
and  measured  in  the  mechanical  stress  experiment 


EXPERIMENTAL  FREQUENCY-TEMPERATURE  CURVES 


b  -  for  the  corrected'cut  ((J)  =  0, 0  =  41.8“,  ijr  =  48.2“)  close  to  the  STC  cut 


4^=0 


•  TEMPERATURE  COEFHCIENT  (units  :ppm.K"') 


Fig.  11 

Theoretical  contour  line  maps  of  temperature  and  planar  isotropic  stress  coefficients 
Thess  curves  are  used  as  a  guide  line  to  experiment  real  quartz  cuts 
in  the  vicinity  of  the  intersection  of  both  zero-stress  and  zero-temperature  effects 

cut  angles :  $  =  0, 37“  <  0  <  47“ 
propagation  direction :  42“  <  ip  <  50“ 

■  =  real  cuts,  experimentally  found  temperature-compensated 
and  measured  in  the  mechanical  stress  experiment 
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In  the  neighboL  ’*’'ood  of  the  P45  cut,  the  sensitivity  of  temperature-compensation  to 
the  propagation  angle  is  about  -3.5  ppm/“C  for  an  increment  of  + 1®  on  ly.  This, 
figure  is  higher  than  for  the  P30  cut  (4>  =  0,  0  =  -10®,  \y  =  30®)  :  -2.5  ppm/®C  for 
Aiy  =  +1°. 

Second-order  temperature  coefficient  is  very  similar  for  the  two  cuts,  i.e.  about 
(-34x10-9  ppm.K-2).  The  theoretical  value  predicted  by  Sinha  [20]  for  the  STC  cut  is 
(-9x10-9  ppm.K-2)  which  does  not  correspond  to  our  experimental  results. 

In  the  neighbourhood  of  this  cut,  the  sensitivity  of  temperature-compensation  to  the 
propagation  angle  \y  is  about  -3.5  ppm/°C  for  an  increment  of  + 1®  on  'P.  This  figure 
is  higher  than  for  the  cut  (cj)  =  0, 0  =  -10®,  ly  =  30®) :  -2.5  ppmy®C  for  Aiy  =  + 1®. 


Force-frequency  effect 

Figures  12,  13  and  14  show  a  comparison  between  experimental  data  and  the 
theoretical  curves  computed  from  the  following  values  (Table  I). 

Stress  compensation  effects  are  clearly  shown  when  the  circular  delay  lines  are 
submitted  to  a  pair  of  orthogonal  forces  of  the  same  magnitude.  But  in  the  case  of  a 
single  diametral  compression,  the  experimental  points  do  not  exactly  match  the 
theoretical  sinusoidal  curve  (Fig.  12).  The  reasons  of  this  discrepancy  might  be 
either  a  systematic  error  due  to  parasitic  stresses  induced  by  the  experimental 
setup,  or  due  to  the  fact  that  the  dimensions  of  the  delay  (9.5  X  4.5  mm2)  are  not 
small  with  respect  to  the  diameter  of  the  plate  (22  mm). 

To  investigate  the  possible  origin  of  the  discrepancy,  the  following  attempts  were 
made : 

-  using  stress  calculation  taking  anisotropy  into  account 

-  increasing  plate  thickness  to  check  the  influence  of  parasitic  bending  forces. 

-  reducing  bias  introduced  by  systematic  calculation. 


21 


F-97.5  MHz 
AF  (Hz/N) 


120 

100 

50 

30 

40 

20 

0 

-20 

-40 

-30 

-30 

-100 

-120 


-140 


$=0  9—2.5  4'«32.5 


s  * 
\ 


LJia'l 


ooo- 


r-97.6  MHz 

AF  (Hz/N) 

* 

120 1 
lOOf 

sol 
60  V 
40l  ^ 
204  ^ 


<t«0  9—5 


i3 


S 


©■C® - - 


-20 
-40 1 
-SOr 

-sol 
-1001 
-120  4- 
-140- 


"T'v  OOO/  \ 

®  w  -©S. - a - -  ?ra-  -ss _ 0-0-  I? 

Of  4  ^Oo  . 


-20  0  20  40  so  SO  100  120  140  ISO  ISO 


-20  0  20  40  SO  SO  100  1  20  140  160  180 


F-98  MHz 
AF  (Hz/N) 


<{i=»0  9=— 7  'P=31 


120 
100 
so 

SOf*^  46 
40l  'v46 


## 
## , 

-20  i  0  ^  /  o 


**** 


i 


-20 
-40 

-so 

-80 

-100 

-120 

-140 


4 


46/ 


O  Q- 


-20  0  20  40  SO  so  100  120  140  ISO  180  ^ 


F-98.4  MHz 

AF  (Hz/N) 

* 

1 

120  6 
loo¬ 
se  t 
606 

406  ' 

20 
0 

-20 
-40 
-60 
-SO 
-100 
-120 
-140 


$=0  9  —  10  ’4'«30 


^^464S6., 


46 


46^ 


46 


'  46 


- 


4S6„ 


-20  0  20  40  SO  SO  100  120  140  160  180 


Fig. 12 

Experimental  force-frequency  effect  on  temperature-compensated  circular  quartz 
delay  lines  with  a  low  sensitivity  to  planar  isotropic  stresses : 

<I)  =  0  -10°  <  0  <  -2.5°  30°  <  ip  <  33° 

single  diametral  compression ;  *  represents  experimental  values 

solid  sinusoidal  curve :  theoretical  values 

.  double  diametral  compression :  o  represents  experimental  values 

dashed  horizontal  line :  theoretical  values 
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Fig.  13 

Experimental  results  of  the  force-freqyency  effect 

solid  lines :  anisotropic  stress  model 
broken  lines :  isotropic  stress  model 
a :  (ST,X)  cut 

b :  («D  =  0,  0  =  41.8“,  ip  =  48.2“)  close  to  STC  cut 
c ;  (4)  =  0,  0  =  -10“,  ip  =  30“) 
a,  b  and  c :  plate  thickness  Imm 
d :  (<E>  =  0,  0  =  -13“,  ip  =  29“)  plate  thickness :  2  mm 
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Fig.  14 

Experimental  results  of  the  force-frequency  effect  in  the  vicinity  ofP45  cut 

solid  lines :  single  diametral  compression 
broken  lines :  double  diametral  compression 
a :  (<I>  =  0,  6  =  39,  ip  =  46)  plate  thickness :  1  mm 

b :  ($  =  0,  0  =  45,  ip  =  47) 

c  :(<!>  =  0,  6  =  39,  ip  =  46)  plate  thickness :  2  mm 

d :  (4>  =  0,  6  =  33,  ip  =  44.4) 
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c)  Improvements  of  the  model  ;  anisotropic  stress  calculation,  effects  of 
plate  thickness 


Introduction  of  anisotropy  in  the  determination  of  the  mechanical  bias 

The  force-frequency  effect  in  a  circular  disk  according  to  the  isotropic  stress 
calculation  is  given  by 

2P  I  .  .  .  .  \  (Oi.'i 

—  *^33^  Q23)cos2p-2  ajjSm2pj 

O 

where  P  is  the  applied  force  along  a  diameter 

d  is  the  diameter  and  t  the  thickness  of  the  disk 

sqjp  saj2  are  the  in-plane  coefficients  expressed  in  the  axes  (aj.ag) 
parallel  to  the  SAW 

P  is  the  azimuthal  angle  between  the  SAW  and  the  force  {Fig.  7). 


The  simple  geometry  of  the  circular  disk  allows  an  exact  evaluation  of  the  quasi¬ 
static  mechanical  bias  at  the  center,  taking  into  account  quartz  anisotropy. 

According  to  Janiaud's  anisotropic  calculation  [16],  the  stress  distribut'on  is  still  a 
simple  radial  distribution  and  its  expression  at  the  center  of  the  plate  is : 


T“,(p  +  40 


p8  /  4A  A  +  B  \ 

ndt  '8-»-e  +  ^  5  ~  +  y  ' 


(25) 


r^(p  +  40 


p8  /  A+B  \ 
ndt  ^  b  —  i,  +  Y  ' 


(26) 


'po  _  Q  (27) 

13  "  ” 

where  p  is  the  azimuthal  angle  between  the  SAW  and  the  force  (Fig.  7) 

W  is  the  propagation  angle  of  the  SAW 

(e,  Y.  8  are  the  following  combinations  of  compliance  constants  expressed 
in  the  axes  (ajO,  a.,o) 


-e  =  4(S,,-S33) 

(28) 

<;  =  Sjj  +  Sjj  -  2Sj2  -  Sgg 

(29) 

V  =  '‘‘5,5  +  533) 

(30) 

8  =  25,3  +  553  +  3533  +  35,, 

(31) 
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Coefficients  A  and  B  are  determined  numerically  by  the  static  equilibrium 
conditions  for  the  plate.  These  coefficients  depend  on  crystal  anisotropy  and  are 
computed  for  each  set  of  cut  angles  and  propagation  direction  T'. 

Final  expression  ofijthe  relative  frequency- shifts  for  a  single  diametral  compression 
of  magnitude  P  is 

y  ~  \  2 

0 

(p+qq-T"  (p+v) 

'“33^  (  - 2 - ) 


*11  +  “33’ 


T°j(p+qq-T^(p-mq 

2 


(32) 


In  that  case,  the  "mean  value"  corresponding  to  a  4-point  symmetrical  compression 
of  magnitude  P/2  is  no  longer  independent  of  the  angles  JP  and  and  is  given  by 
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Comparison  between  isotropic  and  anisotropic  stress  calculation  is  shown  on  Fig.  13 
where  solid  lines  represent  theoretical  values  of  the  anisotropic  stress  calculation 
(broken  lines  —  isotropic  calculation). 

The  difference  between  the  two  models  depends  on  the  crystal  cut.  In  the  case  of  the 
ST-cut  (Fig.  13a)  the  anisotropic  model  predicts  a  well-marked  oscillation  in  the 
case  of  a  4-point  compression.  The  average  value  of  the  oscillation  corresponds  to  a 

constant  value  given  by  the  isotropic  stress  calculation. 

\ 

The  difference  between  the  two  models  is  also  significative  in  the  case,  of  Fig.  13b 
(^  =  0,  0  =  41.8,  W  =  48.2).  In  the  case  of  the  cuts  (4>  =  0,  0  —  -10°,  'P  —  30°)  the 
two  models  predict  almost  the  same  result,  and  do  not  explain  the  discrepancy  in 
amplitude  and  phase  between  theory  and  experiment  on  Fig.  13c.  Fig.  13d  shows  the 
investigation  on  a  thicker  plate  (2  mm  thick)  with  a  better  agreement  between 
theory  and  experiment  ;  this  could  be  explained  by  a  significant  reduction  of 
parasitic  bending  forces  with  respect  to  the  previous  experiments  (1  mm  thick). 


27 


Table  2  summarizes  the  experimental  results  obtained. 
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II  II 
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X  4=0* 

coupling  factor  AV/V 

■Bil 

4.5x10  ^ 

5.8x1 0“"’’ 

THEORY 

beam  steering  angle 

-5.7* 

~  6* 

0 

2-nd  order 
temperature  coefficient 

-34x10“®  K“® 

EXPERIMENT 

sensitivity  to  planar 
isotropic  stresses 

<3x10  .N  ' 

<5x10  .N  ^ 

Table  2 

I 

Through  this  model,  the  last  step,  at  the  present  time,  of  analytical  calculation  has 
been  reached  for  diametral  in-plane  compression. 

As  announced  in  the  second  interim  progress  report  [26]  numerical  analysis  by 
finite  element  method  were  implemented. 


d)  Comparison  with  results  obtained  by  finite  element  analysis 

Up  to  confirm  previous  predictions,  particularly  when  anisotropy  is  held  into 
accoung,  finite  element  model  have  been  coupled  with  the  SAW  calculation 
programm  in  two  manners ; 

-  with  stress  sensitivity  coefficients 

-  with  complete  integral  expression  of  the  perturbation  method. 

Since  the  compression  is  supposed  to  be,  homogeneous  along  the  thickness,  the 
present  model  is  in  two  dimensions  (Fig.  15)'.  As  in  experiments  (Fig.  7),  the  bottom 
of  the  plate  is  fixed  and  forces  are  applied  at  its  top.  The  elastic  constants  used  for 
calculation  are  C,,,  C,y,  C,5,  (D^g,  Ggj,  Cgg  and  the  direction  of  propagation  'F  is 
incremented  with  a  5°  step  up  to  simulate  the  azimuthal  angle  p  between  W  and 
force  direction  as  in  experiments. 
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Fig.  15  :  Two  dimension  finite  element  model 
One  Newton  force  is  applied  at  the  top  of  the  plate, 
the  bottom  remains  fixed  as  in  experiments 


Stresses,  strains  and  displacement  gradients  are  obtained  near  the  center  of  the 
plate,  for  8  triangular  elements,  forming  a  regular  octagon  around  it.  The  value  of 
stress,  strain  and  displacement  gradients  are  known  in  the  barycentre  of  each 
element.  The  degree  of  interpolation  polynomials  is  1  (easier  for  calculating 
displacement  gradients).  After  calculation,  the  datas  have  to  be  rotated  to  bring 
them  back  in  the  propagation  axis. 

In  the  first  case  of  coupling  finite  element  datas  and  the  SAW  programm,  eq.  (18)  is 
used  without  more  complications  if  the  elements  have  the  same  surface.  For  a 
complete  calculation  (2nd  case),  perturbation  method  integral  is  replaced  by 
summation  over  the  elements,  taking  their  surface  into  account.  The  exponential 
term  in  a2  (thickness  coordinate)  is  analytically  integrated  from  0  to  infinite.  In  the 
case  of  the  regular  octagon,  the  8  elements  present  the  same  surface,  so  it  simplifies 
the  expression  used,  but  the  principle  remains  the  same,  as  follows 
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where  e  represents  the  number  of  the  current  element,  S(e)  its  surface  and  Hikjm  ie) 
the  mechanical  prestrained  tensor  for  the  element.  In  the  case  of  isosurface  for 
elements,  the  term  S(e)  is  equal  to  unity. 


Computations  were  made  for  (ST,X)  cut  (Fig.  16). 
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Fig.  16 

Finite  elements  analysis  results  for  (ST,X)  cut 
a :  comparison  with  experimental  measurements 
b :  comparison  with  analytical  isotropic  model 

d)  complete  coupling  between  finite  element  and  perturbation 
methods 

®  calculation  using  stress  sensitivity  coefficients 
@  classical  isotropic  model 


As  in  anisotropic  analytic  model,  a  strong  oscillation  appears  for  (ST,X)  cut  with  the 
complete  reckoning,  that  seems  to  valid  previous  computations.  This  efTect  is  also 
present  when  using  stress  sensitivity  coefficients,  but  very  smoother  than  the 
previous  one.  It  can  be  also  noted  that  the  magnitude  in  the  first  case  is  much 
important  than  experimental  results.  This  might  be  explained  by  strict  restrictions 
on  displacement  used  to  resolve  the  static  mechanical  problem,  and  also  the  fact 
that  interpolation  polynomials  are  of  degree  1,  that  yields  to  significative  errors  for 
the  ihodel.  Presently,  calculation  with  polynomials  of  2nd  degree  is  implemented. 
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As  a  conclusion  of  this  chapter,  one  can  establish  that  analytical  and  numerical 
results  are  in  good  agreement,  and  also  with  most  of  experimental  ones. 

Classical  isotropic  model  can  be  used  to  choose  a  cut  for  its  mechanical  properties. 
More  complete  reckoning  permits  to  improve  simple  models  for  a  given  structure. 
The  isotropic  model  has  been  numerically  implemented  to  test  the  compatibility 
between  finite  elements  datas  and  the  SAW  program.  It  results  that  such  a  manner 
to  calculate  relative  frequency  shifts  seems  better  than  the  stress  coefficients 
calculation  method.  This  phenomena  is  perhaps  due  to  cumulative  errors  along  the 
systematic  tensorial  calculation. 
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4  -  EXPERIMENTAL  STUDY  OF  DYNAMICAL  TEMPERATURE  EFFECTS 


In  this  part,  the  frequency  instability  due  to  dynamical  temperature  perturbations 
is  studied  in  reference  to  previous  works  about  this  subject  [23-24]. 

It  has  been  shown  that  mounting  design  has  an  important  influence  on  the  thermal 
sensitivity.  A  strong  analytical  effort  has  been  made  for  demonstrate  the  existence 
of  cuts  insensitive  to  dynamic  thermal  effects,  and  yields  to  the  conclusion  that 
radial  thermal  conductivity  may  induce  more  frequency  shifts  than  thickness 
conductivity  in  identical  heating  conditions  [24].  Particularly,  a  two-dimension 
isotropic  model  was  established,  giving  the  time  dependent  distributipn  in  a 
circular  plate.  Nevertheless,  the  solution  obtained  was  really  cumbersome  and  not 
of  practical  interest  in  SAW  device  modellings. 


So,  an  experimental  investigation  has  been  presently  choosen  for  the  previous  cuts 
mechanically  tested  (§  3),  with  the  idea  that  dynamic  thermal  perturbation  within 
the  plate  has  a  similar  behaviour  as  for  stress  sensitivity  previously  studied,  but 
with  a  better  accuracy. 


The  SAW  frequency-temperature  variation  law,  given  by  [23] 

^f  ^  dT 

—  =  a  (T~T)+biT-Tr  +  a  — 

f  0  0  0  o  dt 
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where  Uq.  are  respectively  the  first  and  second  order  static  temperature 
coefficients  and  a  the  dynamic  temperature  coefficient,  reduces  to 
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(36) 


when  the  device  is  just  around  temperature  Tq. 

This  model  shows  that  the  relative  frequency  shift  due  to  dynamic  temperature 
effect  is  proportional  to  the  thermal  excitation  speed. 


a)  Experimental  method 

In  fact,  when  a  SAW  resonator  is  submitted  to  dynamic  temperature  biases,  the 
whole  enclosure  is  excited  and  the  thermal  wave  is  assumed  to  propagate  mainly 
via  the  mounting  of  the  quartz  plate.  To  obtain  such  experimental  effects,  a  supply 
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has  been  built  composed  of  a  copper  plate  (choosen  for  its  good  thermal 
conductivity)  with  a  ceramic  case  NPN  transistor  fixed  on  its  bottom  (Fig.  17). 

Two  mountings,  supporting  the  quartz  plate,  are  attached  on  the  top  of  the  copper 
plate  permitting  thermal  excitation  of  the  SAW  propagating  substrate.  The  whole 
system  is  placed  in  a  brass  case  electrically  refered  to  zero  potential  (Fig.  18)  with 
no  thermal  conductivity  between  the  enclosure  and  the  SAW  device. 

Two  resistive  sensors,  one  fixed  under  the  quartz  plate  and  the  other  one  on  one  of 
the  mountings,  give  the  temperature  variations.  This  device  and  the  loop  amplifier 
are  at  the  initial  temperature  of  -30°C  up  to  permit  thermal  cycles  approximatively 
around  ambiant  temperature  (18  to  25°C).  The  experiment  is  controlled  by  a 
programm  that  commands  the  electric  excitation  of  the  transistor  and  stores 
temperature  and  frequency  datas  during  cycles.  Each  experiment  is  made  with  a 
significative  number  of  temperature  cycles,  that  yields  to  obtain  repetitive 
measurements  around  a  steady  state. 


b)  Results 


Classical  (ST,  X)  cut  and  a  single-rotated,  close  to  TG,  cut  (4>  =  0°,  9  =  -7°,  ¥  =  31“) 
have  been  tested,  and  comparison  is  made  between  results.  First,  we  verify  that 
experimental  conditions  of  dynamic  temperature  tests  were  present  by  cycling  both 
cuts  at  different  magnitude  of  thermal  excitation,  for  a  fixed  cycling  time.  Results 
are  shown  in  Table  HI. 

Comparison  between  the  two  cuts  behaviour  is  possible,  and  we  showed  (Fig.  19) 
that  ((4>  =  0“,  6  =  -7“,  V  =31“)  cut  is  about  ten  times  less  sensitive  than  (ST,X)  cut, 
that  seems  to  confirm  previous  assumptions  on  stress  sensitivity. 

It  has  been  demonstrated  that  temperature  behaviour  is  an  accurate  method  to 
determine  the  best  insensitive  cut  for  SAW  device  applications.  Nevertheless,  this 
method  needs  a  heavy  setup  to  obtain  the  sensitivity  and  it  could  be  interesting  to 
develop  numerical  modelings  to  predict  the  sensitivity,  taking  into  account  the 
right  design. 
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Fig.  19  :  Comparison  between  (ST, X)  cut  (a)  and  ((j[)  =  0°,  6  =  -7®, 'F  =  3r)cut(b) 

cycle  duration :  180  s 
number  of  measurements :  171 
Fixation  points  are  the  in-plane  stress  insensitive  ones 


(ST,  X)  cut 

(4,  =  O®;0  =  -7®;'F  =  3r)cut 

cycling  speed 
in  ®K/mn 

(P*f^o)max  • 

in  Hz 

cycling  speed 
in  °K/mn 
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o)min 
in  Hz 

0.203 

19.5 

0.388 

4.3 

0.399 

39.3 

0.433. 

4.6 

0.612 

58.5 

0.788 

9.3 

Tables 


After  the  results  of  this  table,  dynamic  thermal  coefficient  a  of  the  two  cuts  can  be 
evaluated. 

For  ST-cut  a  is  equal  to  1.10  °  s/°K  and  the  (4)  =  0°,  6  =  -7°,  =  31“)  cut  has  a 

dynamic  thermal  coefficient  a  equal  to  .ll.lO'®,  then  ten  times  lower  than  for  the 
usual  quartz  cut. 

This  confirms  the  stress  sensitivity  calculated  in  the  third  chapter. 
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5>  CALCULATIONS  ON  FLEXURAL  SENSITIVITY  OF  SAW  QUARTZ 
CUTS  COMPENSATED  FOR  PLANAR  ISOTROPIC  STKKSSKS 


In  this  section,  the  behaviour  of  a  SAW  device  built  on  a  quartz  plate  submitted  to 
bending  forces  will  be  examined.  The  simple  case  of  a  rectangular  plate,  either 
clamped  at  one  end  and  loaded  at  the  other  (cantilever  beam)  or  clamped  on  two 
edges  and  submitted  to  a  normal  acceleration,  and  the  case  of  the  same  plate 
subjected  to  equal  and  opposite  bending  moments  at  its  ends  are  presented.  In  the 
two  first  approaches,  stress  distributions  will  be  calculated  in  the  classical  model  of 
the  isotropic  thin  plate  [21],  and  their  influence  on  SAW  propagation  will  be 
calculated  according  to  the  model  of  stress  sensitivity  coefficients. 

The  last  one  consists  in  a  complete  anisotropic  modelising  for  stresses,  strains  and 
displacement  gradients,  and  the  numerical  processing  established  to  couple  finite 
element  method  and  SAW  program  will  be  used  to  calculate  the  frequency  shift. 

In  all  these  cases,  the  fixation  points  are  supposed  to  be  "not  too  close"  to  the  SAW 
propagating  beam,  so  boundary  effects  can  be  neglected  (St  Venant  principle). 


a)  classical  isotropic  stress  calculations 


Let  us  consider  the  case  of  the  isotropic  cantilever  beam.  The  width  of  the  plate  is 
supposed  infinite  and  the  loading  force  homogeneous.  So,  in  agreement  with 
classical  thin  plate  theory  [21],  the  problem  reduces  to  the  resolution  of  a  two- 
dimension  differential  equation  which  expression  is 
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with  4)  the  stress  function  related  to  stresses  by  the  relations 
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where  the  absence  of  body  forces  is  assumed  (Fig.  20.). 
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Fig.  20  :  Stresses  in  a  cantilever  beam  clamped  at  one  end 


In  such  a  case,  solution  of  eq.  (37)  in  the  form  of  polynomials  is  of  interest.  Having 
regard  to  boundary  conditions  which  are : 

-  no  stresses  on  longitudinal  sides  ±  h 

-  the  sum  of  the  shearing  forces  distributed  over  the  loaded  end  must  be  equal  to  F 

-Til  must  be  equal  to  zero  at  the  loaded  end  (no  compressive  stress  in  the  ai 
direction  at  the  loaded  end)  and  taking  the  stress  function  in  the  form  of  a 
polynomial  of  fourth  degree  yield  to  the  following  solution 
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By  integrating  these  expressions  over  the  SAW  propagating  volume,  and  then 
dividing  the  results  by  the  value  of  this  volume,  the  mean  values  of  stresses  that 
perturbate  SAW  propagation  are  obtained.  Using  a  perturbation  method  with  the 
hypothesis  of  no  rigid  rotation  of  the  plate  (after  application  of  static  bias),  the 
frequency  shift  is  given  by 
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(40) 


where  <  >  stands  for  the  mean  value  over  the  volume  where  the  SAW  propagates. 


In  eq.  (40)  i;he  loading  force  is  applied  along  the  propagation  direction. 
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The  model,  of  cantilever  beam  has  been  used  for  SAW  accelerometer  applications, 
but  the  case  of  a  rectangular  plate  clamped  on  two  edges  and  bent  by  its  own  weight 
(which  corresponds  to  the  case  of  quasi-static  normal  accelerations)  is  probably 
more  realistic  for  the  modelization  of  SAW  oscillators. 

The  calculation  process  is  the  same  as  the  previous  one,  except  that  a  polynomial  of 
the  fifth  degree  is  chosen  for  the  stress  function  and  the  boundary  conditions  are : 

-  the  plate  is  sunk  at  its  two  ends 

-  there  is  no  normal  stresses  acting  on  the  longitudinal  sides 

-  there  is  no  flexural  moments  applied  at  ±  ^  (Fig.  21). 
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Solutions  respecting  boundary  conditions  given  by  the  classical  theory  of  elasticity 
[17]  yield  the  following  expressions  of  stresses 


(43) 
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The  relative  velocity  shifts  of  the  SAW  will  be  expressed  (using  a  similar  notation 
as  in  eq.  (40)  by 

+’“22  («> 

0 

In  eq.  (44)  the  contribution  of  the  shear  component  <T\2^>  vanishes  by  symmetry. 
Eqs.  (40)  and  (44)  therefore  exhibit  a  similar  form  combining  an  in-plane 
compressive  stress  <Tu°>  and  a  second  term  either  due  to  a  shear  component 
(plate  clampled  at  one  end)  or  to  a  normal  compressive  stress  (plate  submitted  to  a 
normal  acceleration). 

An  important  fact  to  point  out  is  that  if  the  main  contribution  is  due  to  the 
terme  squ  <Tii>,  flexural  sensitivity  will  be  equivalent  to  in-plane  stress 
sensitivity.  Under  this  assumption,  a  device  built  with  a  cut  compensated  for 
planar  isotropic  stresses  will  also  exhibit  a  low  sensitivity  to  a  symmetrical 
distribution  of  flexural  stresses. 


b)  Anisotropic  stress  model 


Timoshenko  and  Goodier  [21]  established  that  a  3-dimensional  elasticity  problem 
can  be  modelised  with  equilibrium  equations  and  boundary  conditions  for  stress 
components  either  constant  or  linearly  dependent  of  coordinates,  so  the  compati¬ 
bility  equations  are  identically  respected,  and  the  stresses  are  correct  solutions  of 
the  problem.  If  the  stress  distribution  is  assumed  to  be  as  follows 


(45) 


with  I  the  moment  of  inertia,  M  the  bending  moment  and  02  the  thickness 
coordinate  (Fig.  22),  we  showed  that  anisotropy  could  be  introduced  for  calculating 
strains  and  displacements  without  complication  for  the  model. 
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The  reckoning  procedure  is  quite  the  same  as  exposed  in  [21]  and  yields  to  the 
following  results : 
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where  conditions  of  fixed  point  in  the  center  of  the  plate  and  no  rigid  rotations 
around  it  have  been  imposed  up  to  assure  the  unicity  of  the  solution. 


Values  for  stresses,  strains  and  displacement  gradients  have  been  calculated  in  the 
propagation  area  for  a  certain  number  of  small  parallelepipedic  elements  (Fig.  23), 
and  then  introduced  in  the  perturbation  method.  Fig.  26  shows  the  results  for 
(ST,X),  experimental  STC  and  TG  (-10,  30)  cuts,  which  seems  to  confirm  the 
possibility  of  very  low  sensitivity  to  flexural  effects  and  especially  for  TG  cut. 
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c)  Numerical  results  of  the  cuts  (P30),  (P45)  and  (ST«X) 

The  flexural  effect  was  calculated  for  a  device  built  on  a  rectangular  plate  with 
some  realistic  values 

plate  lenght :  17  mm  -+  22  mm 
thickness :  2  mm 
wavelength :  34.4  pm 
wavelength/thickness  ratio :  1.7  X  10'^ 

Results  are  presented  on  Fig.  24  for  the  plate  clamped  at  one  edge  and  on  Fig.  25  for 
the  plate  clamped  on  two  edges  and  submitted  to  a  normal  acceleration.  Results  are 
expressed  in  terms  of  relative  velocity  shifts  (normalized  to  the  mean  value  of 
compressive  strains  <Sii>  =siin<Tii>)asa  function  of  the  azimuthal  angle  p 
between  the  SAW  propagation  direction  and  the  edges  of  the  rectangular  plate 
(Figs.  24a  and  25a). 

For  the  P30  and  P45  cuts,  the  main  contribution  to  the  flexural  effect  is  due  to  the 
compressive  stress  <T\\>,  The  sinusoidal  shape  of  the  curves  on  Figs.  24  and  25 
may  be  interpreted  siniply  as  the  tensorial  transformation  of  stresses  when  turning 
the  azimuthal  angle  p  ;  this  behavior  is  very  similar  to  the  case  of  simple  diametral 
compression  calculated  in  sec.  3,  Eq.  (22). 

Results  obtained  with  the  anisotropic  model  appear  on  Fig.  26.  They  show  a  strong 
similarity  about  the  shape  with  the  anisotropic  radial  in-planes  sensitivity  curves 
(section  3)  and  a  general  identical  behaviour  with  previous  isotropic  models 
(because  of  the  particular  stress  distribution,  that  is  T\\  the  only  non-zero  stress 
term). 
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Fig.  24 

Flexural  effects  (isotropic  stress  model)  for  a  SAW  propagating 
on  a  bent  cantilever  beam 

dimensions :  L  =  17  mm  thickness  =  2  mm  K  =  34.4  pm 
a :  definition  of  the  azimuthal  angle  P 
b,  c,  d :  comparison  of  planar  isotropic  stress  compensated  cuts 
with  the  classical  (ST.X)  cut 
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Fig.  26 

Flexural  effects  for  a  SAW  propagating  on  a  quartz  plate 
submitted  to  equal  and  opposite  bending  moments  M. 

These  results  are  very  closed  to  previous  bending  effect  models  and  to  compressive 
in-plane  anisotropic  model  because  of  the  particular  stress  distribution 
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For  a  suitable  choice  of  the  angle  P,  the  device  will  be  insensitive  to  flexural  effects. 
In  the  case  of  the  P30  and  P45  cuts,  there  is  no  offset  value  of  the  sinusoidal  curve, 
and  therefore  the  device  will  be  insensitive  to  a  symmetrical  distribution  of  bending 
stresses  :  this  will  happen  if  the  two  principal  compressive  stresses  T\i  and  Tzz  are 
equal  and  with  same  sign  ;  the  contribution  of  other  stresses  such  as  T12,  T32,  T22 
being  of  much  lower  magnitude. 


CONCLUSION 


Experimental  measurements  have  been  extended  in  the  vicinity  of  the  P45  cut, 
close  to  B.K.  Sinha's  STC-cut  [14].  The  following  cuts  ((|)  =  0“,  6  =  39°,  'P  =  46°), 
((])  =  0°,  B  =  33°,  'P  =  44.4°),  (4)  =  0°,  6  =  45°,  'P  =  47°)  have  been  experimentally 
found  to  be  temperature-compensated  (by  small  correction  on  the  propagation 
direction  'P  relatively  to  the  theoretical  values  given  by  computation)  and  then 
mechanically  tested  in  radial  in-plane  compression  experiment.  During  these 
manipulations,  it  has  been  shown  as  previously  than  reducing  stress  sensitivity  by 
increasing  thickness  of  the  plate  induces  higher  experimental  dispersion  of  the 
results. 

Finite  elem.ent  method  has  been  employed  in  the  calculation  process  of  frequency 
shift.  It  was  of  major  importance  to  valid  the  previous  analytical  anisotropic  model. 
Results  that  have  been  obtained  are  rather  in  a  good  agreement  with  experiments 
for  a  two  dimensions  model.  More,  the  tests  implemented  in  order  to  establish  the 
validity  of  the  coupling  between  finite  elements  and  SAW  program  showed  that 
systematic  calculations  may  introduce  significative  errors  in  the  stress  sensitivity 
coefficients  calculation. 

A  dynamic  thermal  experiment  was  implemented.  It  has  been  shown  that, 
depending  on  SAW  orientations,  it  exists  particular  cuts  which  are  less  sensitive 
than  others,  or  particularly  than  standard  (ST,X)  cut.  That  may  yield  to  the 
realisation  of  SAW  devices,  insensitive  to  in-plane  compression,  flexion  effects, 
statically  and  dynamically  temperature-compensated,  with  a  satisfying  safety 
margin  on  cut  angles  especially  in  the  vicinity  of  the  P30-cut. 

A  new  flexural  anisotropic  model  was  established  with  a  good  agreement  between 
its  results  and  the  isotropic  previous  one,  given  in  the  2nd  interim  progress.  More, 
since  St  Venant's  principle  is  assumed  so  that  only  compressive  effects  are  of 
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importance  in  the  plate  center,  the  accordance  between  all  the  mechanical  models  is 
clearly  demonstrated. 

Of  course,  we  must  purchase  our  investigations,  particularly  for  dynamic  thermal 
effect  experiments  and  computations,  that  are  now  possible  by  the  use  of  finite 
element  program  MODULEF  (by  courtesy  of  Mr.  Crolet,  Director  of  the  L.C.S.; 
besangon,  France).  We  will  investigate  the  time  dependent  thermal  distribution  in 
a  3  dimensions  anisotropic  plate  and  the  stresses  due  to  a  dynamic  temperature 
excitation. 

We  will  also  attempt  to  work  on  higher  precision  quartz  devices,  up  to  confirm  the 
analytical  predictions,  and  cuts  in  the  vicinity  of  P45  cut  will  be  thermally 
experimented. 

Also  classical  quartz  lines  used  for  the  previous  experiments  should  be  replaced  by 
real  resonators,  that  would  introduce  some  modifications  in  the  theoretical 
approach  of  phenomena. 

The  different  points  can  merge  into  a  general  method  to  define  the  whole  structure 
of  a  low  sensitive  SAW  device  for  military  or  space  applications. 
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APPENDIX  1 

CUT  ANGLES  AND  PROPAGATION  DIRECTION  FOR  SAWs  ON  QUARTZ 
DEFINITION  AND  SYMMETRY  PROPERTIES 


Fig.  27  ;  Correspondance  between  IEEE  and  Euler  angles 


Correspondance  between  the  two  angles  is  obtained  by  identifying  the  unit  vectors 
(ai,  a2,  as)  defined  by  IEEE  angles  and.(B,  N,  F)  vectors  using  Rosenberg's  notation 
for  Euler  angles. 
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The  only  difference  appears  in  the  definition  of  the  second  angle  0ieee  or  p  Euler,  the 
first  refering  to  the  plate  surface  and  the  second  to  the  plate  normal.  Thus,  the 
correspondence  is  simply : 

^Euler  =  ‘^ieee 

PEuler  ~  ®iEEE  *  (A.2) 

0Euler  =  '^^ieee  . 

2  -  Symmetry  properties  (IEEE  notation) 

Surface  waves  on  quartz  obey  the  three  following  symmetries. 

a)  Redundancy  in  the  definition  of  the  entangles 

From  the  definition  of  the  cut  angles,  the  following  symmetry  property  appears : 

(4*.  0,  ip)  =  (4)  +  180°,  180°-  0,  'F  +  180°)  (A.4) 

b)  Symmetries  of  quartz 

•  3-fold  symmetry  around  the  Z-axis 

(4),  0,  ip)  =  (4)  +  120°,  0, 'F)  (A.5) 

•  2-fold  symmetry  around  X-axis 

(4),0,ip)  =  (-4),0  +  180°,W)  (A.6) 

Combination  of  (29)  and  (27)  yields  an  equivalent  form 

(4),  0,  ip)  "  (180°-4),  -0,  180°  +  W)  (A.7) 

c)  Symmetries  of  surface  waves 

Two  surface  waves  on  the  same  plate  have  the  same  properties  in  case  of : 

•  propagation  on  the  same  surface,  with  opposite  directions 

(<^,e,ip)  =  ($,0,W  +  18O°)  (A.8) 

•  propagation  in  the  same  direction,  but  on  the  other  side  of  the  plate 

(4>,e,ip)  =  (4>,0  +  18O°,-ip) 


(A.9) 
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Eqs.  (A.8)  and  (A.9)  come  from  the  symmetry  of  the  determinants  that  define  the 
surface  wave  parameters. 

Combination  of  (A.6)  and  (A.9)  yields  the  equivalent  form 

(4>,  6,  W)  =  (A.IO) 


d)  Consequences  on  the  angular  range  required  for  a  complete  stress  and 
temperature  coefficients  mapping 

Starting  from  the  most  general  angular  range  :  0  -►  360°  for  the  three  angles,  the 
three  symmetries  a)  b)  and  c)  yield 

-  eq.  (A.4)  reduces  the  angular  range  on  0  to  -90°  -►  +  90° 

-  eq.  (A. 5)  reduces  the  angular  range  on  ^  to  -60°  -♦  +  60° 

-  eq.  (A.7)  reduces  the  angular  range  on  8  from  -90  -»  +  90°  to  0  -♦  +  90° 

-  eqs.  (A.8)  and  (A.9)  reduce  the  angular  range  on  qr  to  0  90°. 

The  combination  of  (A.5),  (A.7),  (A.8)  yields : 

(4),0,v{/)  =  (60°.(i),-e,'P)  (.4.11) 

(4),  0,  \j()  =  (-60°-(j[), -0,  ^)  (A.12) 

Eqs.  (A.ll)  and  (A.12)  show  that  the  points  (({>  =  +30°,  0  =  0)  and  (4)  =  -30°,  0  =  0) 
are  centers  of  symmetry  for  the  contour-lines  on  Figs.  4  and  5. 


3  -  Correspondance  between  IRE-1949  and  IEEE-1978  standards  on 
piezoelectricity 

In  1978,  the  IEEE  standard  on  piezoelectricity  [27]  changed  the  orientation  of  the 
canonical  axes  of  quarti,  so  that  the  definition  for  quartz  should  be  coherent  with 
other  crystals.  The  consequence  of  this  change  is  a  rotation  of  the  canonical  axes  by 
180°  around  the  Z-axis  i.e. :  X  ->  -X,  Y  -♦  -Y. 

Most  of  authors  still  refer  to  the  "old"  1949-standard  which  is  taken  as  the  reference 
in  this  work.  Correspondance  between  the  two  standards  is  as  follows  :  a  cut  defined 
by  the  angles  (4>,  8,  qO  in  the  1949-standard  will  be  the  same  if  defined  by  the  angles 
(4) +  180°,  0,  ip)  in  the  1978-s^ndard.  Using  symmetry  equation  (A.7)  the 
correspondance  between  the  two  standards  is  finally : 
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1949-cut  ((i>,  0,  Ilf)  =  1978-cut  (-(^,-0,iif)  (A.13) 

For  example  :  the  ST-cut  will  be  defined  by  (4)  =  0,  0  =  42.75°,  'P  =  0)  in  1949- 
standard  and  by  (4)  =  0, 0  =  -42.75°,  'F  =  0)  in  the  "new"  1978-standard. 
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Force  Command,  Control,  Communications  and  Intelligence  (C  l)  activities 
for  all  Air  Force  platforms.  It  also  executes  selected  acquisition  programs 
in  several  areas  of  expertise.  Technical  and  engineering  support  within 
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ESD  elements  to  perform  effective  acquisition  of  C  l  systems.  In  addition, 
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Air  Force  user  community,  and  other  DOD  and  non-DOD  agencies.  Rome 
Laboratory  maintains  technical  competence  and  research  programs  in  areas 
including,  but  not  limited  to,  communications,  command  and  control,  battle 
management,  intelligence  information  processing,  computational  sciences 
and  software  producibility,  wide  area  surveillance/sensors,  signal  proces¬ 
sing,  solid  state  sciences,  photonics,  electromagnetic  technology,  super¬ 
conductivity,  and  electronic  reliability/maintainability  and  testability. 


